Bichromatic Slowing of Metastable Helium 
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We examine two approaches for significantly extending the velocity range of the optical bichro- 
matic force (BCF), to make it useful for laser deceleration of atomic and molecular beams. First, 
we present experimental results and calculations for BCF deceleration of metastable helium using 
very large BCF detunings, and discuss the limitations of this approach. We consider in detail the 
constraints, both inherent and practical, that set the usable upper limit of the BCF. We then show 
that a more promising approach is to utilize a BCF profile with a relatively small velocity range 
in conjunction with chirped Doppler shifts, to keep the force resonant with the atoms as they are 
slowed. In an initial experimental test of this chirped BCF method, helium atoms are slowed by 
~ 370 m/s using a BCF profile with a velocity width of < 125 m/s. Straightforward scaling of 
the present results indicates that a decelerator for He* capable of loading a magneto-optical trap 
(MOT) can yield a brightness comparable to a much larger Zeeman slower. 

PACS numbers: 37.10.De,37.20.+j 



I. INTRODUCTION 

Most laser slowing methods rely on the radiative force 
in a monochromatic field [H-Q- The momentum transfer 
Hk arising from photon absorption allows atoms to be ac- 
celerated or decelerated, and their velocity distribution 
to be narrowed (cooled). However, the radiative force is 
limited in its magnitude by the spontaneous emission rate 
7, to a maximum of i^ad = hkj/2. Long slowing times 
allow the small transverse velocity spread in an atomic 
beam to significantly reduce the source brightness, neg- 
atively impacting the loading of the slowed atoms into a 
magneto-optical trap (MOT). 

Unlike the radiative force, the dipole force arising from 
intensity gradients in a monochromatic standing wave is 
not limited by the radiative decay rate. However, it al- 
ternates sign on the scale of an optical standing wave. 
The average force is zero, limiting its usefulness in a de- 
celerator, although a pulsed slower can be created by 
utilizing a transient standing-wave lattice produced by 
an intense pulsed laser Q . Another way around this lim- 
itation, demonstrated by Kazantsev and Krasnov [||, is 
to rectify the dipole force. Adding a second frequency 
to the standing wave results in a modulated light shift, 
which can be adjusted to keep the sign of the dipole force 
positive. This technique eventually resulted in the de- 
flection of sodium atoms by a few m/s by Grimm and 
coworkers 

Continued refinement of slowing techniques in bichro- 
matic fields by Grimm and coworkers eventually led to 
the observation of a much larger rectified dipole force in 
cesium, with a velocity range of 225 m/s [|. This was the 
birth of the optical bichromatic force (BCF), which relies 
on a pair of counterpropagating two-color beams that, 
when the intensity and standing wavelength are carefully 
selected, will coherently drive the atom through cycles of 
photon absorption and stimulated emission much more 
rapidly than the radiative decay rate. 

These developments led to more extensive research into 



the nature and origins of the bichromatic force, which to 
date has been demonstrated in Cs Q, Na 0, Na 2 0, 

Kb [HEl, and He 0~E3- Des P ite this 

progress and the 

obvious advantages of a strong continuous optical force, 
an atomic beam slower for MOT loading based on the 
BCF has never been realized. 

We have designed two experiments to test the poten- 
tial limits of the BCF and fully realize an atomic beam 
slower. The first set of experiments tests the sensitivity 
and limitations of the BCF at large detunings. This is 
required for a static atomic slower design, which would 
utilize one or two pairs of bichromatic beams configured 
such that the BCF profile has a wide velocity range, ca- 
pable of decelerating He* atoms by several hundred m/s. 

The second set of experiments explores an alternative 
BCF slower scheme that avoids many of the difficulties of 
using large bichromatic detunings. Here BCF beams with 
a relatively small detunings are used, with the addition of 
dynamic frequency chirping to maintain resonance with 
the changing Doppler shifts of the decelerating atoms. 
This design conserves laser power at the cost of added 
complexity in the control electronics. More importantly, 
it keeps the BCF parameters in a range within which they 
are robust against small misadjustments. We demon- 
strate deceleration by 370 m/s in a prototype design, and 
discuss straightforward improvements that should double 
this range. 



II. MODELING THE OPTICAL BICHROMATIC 
FORCE 

A. 7r-pulse model 

An intuitive but somewhat inaccurate model for the 
BCF is based on the concept of periodic momentum 
exchange between the atom and the co- and counterprop- 
agating bichromatic beams using 7r-pulses. As shown in 
Fig. 1, each pair of two-color beams is offset by the 
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FIG. 1. Simple 7r-pulse model of the optical trichromatic force 
on an atom, adapted from Fig. 1 of Ref. [l8| . Two-color 
beams impinge from each direction, offset by frequencies ±<5 
from resonance at u a . For large atomic velocities v, Doppler 
offsets of ±fcu must also be added as shown. The frequency 
pairs at ±5 give rise to a series of rf beat notes, each with an 
area of approximately n. 

atomic Doppler shift kv so that in the frame of the mov- 
ing atom, the frequencies are offset from the atomic res- 
onance ujq by ±<5. Thus the atom experiences from each 
direction a traveling sequence of beat notes, each with 
duration T = it/ 5. In the 7r-pulse model, the two sets 
of beat notes are treated as if they are independent non- 
overlapping pulse trains. In this approximation, if the 
laser power is chosen so that the beat note area satisfies 
the 7r-pulse condition for an atomic transition, the atom 
will experience a sequence of alternating excitations and 
stimulated de-excitations. This occurs when 

n = 71-6/4, (l) 

where 51 = 7(//2/ sat ) 1 / 2 is the on-resonance Rabi fre- 
quency of each bichromatic frequency component |2j. A 
large decelerating force to the left occurs if the excitations 
come from the right-hand beam and the de-excitations 
from the left. If 6/n 3> 7, this force will greatly exceed 
the radiative force i^ad- 

This force is quite tolerant of deviations from the ex- 
act 7r-pulse condition so long as the left-hand beat notes 
closely match the right-hand ones. In the pulse-pair 
model, the two pulses will then appear to the atom as 
time-reversed images of one another. The second pulse 
in each pair reverses the curve on the Bloch sphere traced 
out by the first [l9| , leaving the atom back in the ground 
state after each two-pulse cycle. This simple time-domain 
picture is very useful for estimating the breakdown of the 
BCF at large detunings or large atomic velocities, a sub- 
ject that we discuss further in Section lill CI 

The 7r-pulse model also leads to a simple and intuitive 
estimate for the velocity range of the BCF, which appears 
to have been overlooked in the previous literature. To a 
first approximation, finite atomic velocities merely induce 
Doppler shifts +kv and —kv for the two pulses in each 
pair, and since the pulses still appear as time-reversed 
images there is little effect. However, the generalized 
Rabi frequency is also modified , because the Doppler 
shifts detune the center or "carrier" frequency of each 
beat note from the atomic resonance: 

n' = ^n 2 + (kv) 2 . (2) 



If the Doppler shift kv is comparable to the resonance 
Rabi frequency Q, the Rabi cycling even for a single pulse 
is so badly out of phase that the BCF will be disrupted. 
In the simplified limit of a rectangular pulse with O w 5, 
the resulting full width of the velocity range is 

A««2p (3) 

in agreement with much more sophisticated models apart 
from the factor of 2. 

The relative phase between the co- and counterprop- 
agating beat notes is critical to the optimization of the 
BCF, as can easily be seen by considering a single pair 
of 7r-pulses separated by time At, the first arriving from 
the right and the second from the left. If the atom is 
excited by the first pulse but undergoes spontaneous de- 
cay before the second pulse arrives, it will start cycling 
in the wrong sequence, leading to acceleration. This can 
be minimized by making At as short as possible while 
maintaining the discreteness of the pulses. Further, if 
the repetition period T 3> At, an atom that is already 
cycling in the wrong sequence will be more likely to radi- 
ate to the ground state after the two-pulse sequence than 
in between the pulses, placing it back into the correct se- 
quence for deceleration by future pulse pairs. For short 
rectangular pulses the optimal configuration would be 
with minimal delay At, but for beat note trains, the best 
one can do is to set the relative phase to cf) = tt/2, where 
4> = corresponds to the limiting case of simultaneous 
arrival, where the pulses fully overlap and there is no net 
force, and 4> = tt is the equally ineffective symmetric case 
where At = T/2. Averaging over several radiative life- 
times with cj) = 7r/2, the correct (decelerating) sequence 
occurs during 3/4 of the interaction time, while the atom 
is accelerated during the remaining time. Numerical cal- 
culations and experiments support this argument, which 
reduces the average force by a factor of two compared 
with perfect in-phase cycling. Because the individual ir- 
pulses exchange momentum increments of hk at rate 5/tt, 
the average bichromatic force is 

F h = hk5/-K. (4) 

This average force exceeds the radiative force by a factor 
of 2<5/(7T7). 

B. Doubly-dressed atom model 

A more complete model that still provides some intu- 
itive insight involves dressing the atomic levels with the 
red- and blue detuned photon field number states. Orig- 
inally developed by Grimm and coworkers [20L |21| , this 
"doubly dressed" model results in a ladder of states sep- 
arated by hS and an infinite tridiagonal Hamiltonian ma- 
trix that is well described by Metcalf and Yatsenko [22j ■ 
Numerically solving for the eigenstates of the Hamilto- 
nian provides an estimate of the BCF magnitude that 



agrees with the 7r-pulse model, and is also in agreement 
that the optimal phase is ir/4. Unlike the 7r-pulse model 
it predicts that the optimal Rabi frequency is given by 

£V = y/3/2d, (5) 

This result, larger than Eq. [I] by about 56%, is in excel- 
lent agreement both with experiment and more complete 
numerical calculations. Applying Floquet theory and us- 
ing a second-order perturbation treatment, this model 
also yields an estimate of the velocity range of the force, 
but this prediction is too small by about a factor of two, 
probably because the effects of Doppler shifts are not 
fully incorporated flij . 



C. Numerical solutions of the optical Bloch 
equations 

A more comprehensive but less transparent approach 
is to numerically solve the optical Bloch equations 
(OBEs) [H, Hi] in the bichromatic field. The OBEs 
fully describe the internal state of a two-level atom in 
the rotating- wave approximation (RWA), taking into ac- 
count the Rabi cycling of atomic states in the field as well 
as spontaneous emission, which is added phenomenolog- 
ically. We obtain numerical solutions using a program 
based largely on code developed for earlier BCF studies 
by the authors of Refs. [1, [2(J The rotating frame 
is chosen to be at the center frequency of the bichro- 
matic spectrum, leaving some residual time dependence 
due to the approximately symmetric detunings ujq ± 8. 
The OBEs are solved over small time steps in the bichro- 
matic field and the instantaneous force is evaluated using 
Ehrenfest's theorem Q, 

F(z,t) = h[u(t)V$t(n r (z,t)) - v(t)VS(Q r (z,t))], (6) 

and then averaged over several radiative lifetimes. The 
atomic velocity is not explicit in the force calculation, 
but is inferred using a constant-velocity approximation 
so that the atomic position is given by z — vt. This 
process is repeated over a range of atomic velocities to 
produce force profiles like the one shown in Fig. [2] 

At the center of the force profile, the calculated am- 
plitude is maximized when the Rabi frequency satisfies 
Eq. [5] and the bichromatic beat phase offset is <f> = tt/2, 
verifying both of the central predictions of the doubly- 
dressed atom model. Both predictions also agree well 
with experimental results. Under these optimal condi- 
tions the maximum force is -Fbcf ~ fokS/w, also consis- 
tent with both the simplified models and experiments. 

The velocity range of the force in the profile of Fig. [2] 
is about 5 /2k, but the profiles take on distinctly differ- 
ent shapes for 8 < O < ^3/28 [25|], developing broad 
peaks near ±8 /2k with a force nearly as large as for the 
central zero-velocity peak. Atoms traversing a Gaussian 
bichromatic beam profile will interact with a range of 
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FIG. 2. Calculated bichromatic force vs. velocity for 8 = 1547 
and phase (f> = tt/2, for Q pt = \/3/2 8. This optimum Rabi 
frequency yields the largest central force, in agreement with 
Eq. [S] rather than Eq. [T] Note that the force magnitude is 
shown as a multiple of Mry, which is twice the radiative force 

-^rad • 



intensities, so a realistic simulation should include an av- 
erage over the force profiles. This yields an improved 
estimate of the BCF velocity range [m [l7j. 



Av » S/k, 



(7) 



that is half as large as the rough estimate in Eq. [3] 

The time required to slow an atom across this full ve- 
locity range is ieCF = rnAv/FBCF- This "bichromatic 
slowing time" is independent of the detuning 8 and for he- 
lium its value is £bcf = 5.9 /is. ft follows that even a fast 
beam of metastable helium atoms with (v ~ 1000 m/s), 
can be decelerated by Sv in an interaction length of less 
than 1 cm. 

Another feature of the BCF profile is the sharp drop-off 
at either end of its range, clearly visible in Fig. [5J If the 
force is applied for a time greater than the bichromatic 
slowing time, atoms within the velocity range of the force 
will "pile up" at the low- velocity limit. This can yield 
significant cooling in addition to deceleration, as can be 
seen in Fig. H and Refs. H El [13]. Recently Harold 
Metcalf has argued that the cooling is due to radiation 
exchange between the bichromatic fields, and may not 
be subject to limitations tied to the rate of spontaneous 
emission 1261] . 



III. LARGE-DETUNING BCF 

A practical atomic beam slower for metastable helium 
must reduce the atomic velocity by at least 700-800 m/s 
to produce a large flux of atoms at MOT capture veloc- 
ities. Ideally, this could be achieved with a single pair 
of counterpropagating beams producing a BCF profile 
with a correspondingly large velocity range. However, 
this would require unreasonably high irradiances, so it 
is necessary to break up the required velocity range in 
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some way. In this section we investigate the feasibility 
of a static two-stage slower and more generally, we study 
the large-detuning limits of the BCF. 



A. Experiment 

For efficient deceleration by two successive BCF in- 
teractions, the first stage should have a velocity profile 
Doppler shifted so that its leading edge is located in an 
appropriate portion of the slow tail of the atom veloc- 
ity profile, and the second stage should have a smaller 
Doppler shift with its leading edge just overlapping the 
trailing edge of the first profile. We have constructed a 
decelerator in the first-stage configuration to test the con- 
cept, as diagrammed in Fig. [3J The metastable helium 
source is a reverse-pumped, liquid nitrogen cooled dc dis- 
charge source patterned closely on a prior design [13, HH , 
with external modifications to permit installation in an 
existing vacuum chamber. The metastable atom flux 
from the source is about 3 x 10 13 He* atoms/sr-s with 
a most probable velocity of approximately 1050 m/s and 
a metastable fraction of roug hly 5 x 10" 5 [U. The atomic 
beam passes from the source chamber through a 500 /im 
diameter skimmer aperture into a time-of-flight chamber. 
After passing through a 70 ^m diameter collimating aper- 
ture, it is mechanically chopped by a tuning fork chop- 
per operating at 160 Hz with a 100 ^m slit width. The 
resulting chopped atom beam, with a divergence half- 
angle of 4.1 mrad, travels through the counterpropagat- 
ing bichromatic decelerating laser beams for a few cen- 
timeters. After an additional time-of-flight (TOF) path, 
the chopped beam impinges on a stainless steel Faraday 
cup. Electrons ejected by He* atom impacts are collected 
and detected by a Ceratron continuous dynode electron 
multiplier. The He* velocity distribution is then calcu- 
lated from the TOF spectra, with a velocity resolution of 
about 67 m/s. 

The bichromatic laser beam pairs are generated from a 
single diode laser (Toptica DL100) locked to the 2 3 S\ — 
2 3 P2 transition at 1083.3 nm using a saturated absorp- 
tion spectrometer. This laser at the atomic frequency 
uj a is double-passed through an acousto-optic modulator 
(AOM) operating at the bichromatic detuning frequency 
S/2tv, producing two superimposed frequencies at ui a ± S. 
It is then coupled into a cw fiber laser amplifier (Nufern 
NuAmp) with a single-mode output of up to 7 W. The 
amplified beat-note train is split into equal parts, then 
Doppler frequency shifts of ±kv/2ir = ±800 MHz are 
added using additional AOMs, to center the bichromatic 
force profile at a velocity of about 866 m/s. An optical 
delay line is used to set the relative phase between the two 
bichromatic beat note trains to the optimal (f> = 7r/2 
The bichromatic beams enter the chamber linearly po- 
larized, pass through quarter-wave retarders to produce 
tr + light, and cross the helium atomic beam at an angle 
of about 1°. The beams are tightly focused to waists 
with top-hat radii of 0.32 mm to provide the required 



irradiance without exceeding the capabilities of the am- 
plifier. Because of this tight focus the overlap with the 
atomic beam is imperfect, resulting in an interaction re- 
gion about 3.7 cm long but covering only about 40% of 
the atomic beam. 



B. Results and analysis 

The results for this large-detuning slowing are shown 
in Fig. |4] We extend previous work using a fixed detuning 
of 184 7 (-300 MHz) [3 [3 by exploring several differ- 
ent detunings. The figure shows raw experimental data 
without enhancement to compensate for imperfect beam 
overlap, as was done in some prior BCF experiments on 
He* [3 E3 • Consequently, the depth of the holes in the 
velocity profiles can be at most 40% of the population, 
which may be further reduced because the atomic beam 
contains some atoms in the 2 1 Sq metastable state that 
cannot be slowed. 

Unfortunately these experiments were terminated 
abruptly by the failure of our 7 W fiber amplifier. As 
a result, only one set of data was acquired in the previ- 
ously unexplored range of detunings between 1857 (300 
MHz), where a large BCF was observed [l5j], and 3687, 
where no effect could be found 16]. Our new results 
are consistent with these prior reports, showing a badly 
disrupted BCF at an intermediate 2787, establishing an 
appareent upper limit to the BCF that we discuss further 
in Sec. iHICl 

The changes in the velocity profiles shown in Fig. [4] 
were confirmed to be due to the BCF by adjusting the 
beat note phase, by blocking one of the two bichromatic 
beams, and by blocking a single frequency component in 
one of the beams. As expected, when the phase is in- 
creased to 4> — 3n/2, the sign of the force changes and 
the deceleration becomes an acceleration. This is per- 
haps the most obvious and identifying signature of the 
BCF. When the phase is set to = or </> = tt, the force 
vanishes altogether, also as expected. If one or more of 
the frequency components is blocked the BCF vanishes 
and smaller peaks appear in the atomic velocity distri- 
bution due to the radiative force from each bichromatic 
frequency component. 

To establish a consistent convention for measuring the 
extent of slowing and the velocity range of the BCF, we 
define "slowing" as the velocity difference between the 
peak of the slowed atom distribution and the correspond- 
ing minimum of the hole in the slowed velocity distribu- 
tion, as indicated in Fig. The velocity range of the 
BCF can be estimated by measuring the full width be- 
tween the half-maximum of the slow side of the slowed 
atom peak and the half-minimum on the fast side of the 
hole in the original velocity distribution. This is the con- 
vention used by the Metcalf group in their BCF slowing 
experiments (la.[l7|. In some circumstances this may dif- 
fer from the actual full width of the BCF velocity profile, 
but we will adopt this convention to facilitate compar- 
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FIG. 3. (color online). Experimental configuration for observing the BCF in atomic He using detunings of 123 — 2787 with a 
static Doppler shift centered at atomic velocity v. Polarizing beamsplitter cubes are denoted by PBC, and a beamsplitter by 
BS. 
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FIG. 4. (color online). Experimentally observed velocity pro- 
files for BCF slowing, increasing in detuning from top to bot- 
tom. The dotted black trace is the unslowed (background) 
distribution, solid red is the slowed distribution, and solid blue 
is the difference. Notice the anomalous result at 5 — 2787, 
consistent with previous observations of an upper limit to the 
BCF detuning. 



isons between prior results and ours. 

Examining the slowing data summarized in Table U 
we verify that the BCF slowing is proportional to 6, and 
that the velocity range of the force is very close to 8/k 
as predicted by the numerical calculations described in 
Section [Til This agrees wi th p revious results from the 
Metcalf group at 8 = 1857 W- 

Because the length of the slower is about 3.7 cm, the 
atoms experience 6-8 BCF slowing time intervals £bcf, 
depending on the detuning. This complicates a direct 



TABLE I. Summary of BCF slowing results, using the mea- 
surement conventions described in the text. Note that the re- 
sult for S = 2787 does not exhibit a velocity profile consistent 
with unidirectional BCF slowing, so this slowing measurement 
is shown only for comparison. 



Detuning 


Slowin 


y 


Velocity Range 


7 


■y/k 


S/k 


5/k 


123 


67 


0.55 


1.04 


154 


86 


0.56 


1.08 


185 


100 


0.54 


1.00 


278 


46 


0.16 


N/A 



measurement of the BCF magnitude since the atoms will 
be slowed by the full velocity range even when the force is 
well below the optimal value of KkS/2. Attempts to mea- 
sure transient deceleration by pulsing the bichromatic 
beams were not successful because of the limited time 
resolution of the chopped-beam TOF scheme. 

There is much insight to be gained from the anomalous 
velocity profile obtained with the largest BCF detuning 
of 5 = 2787 (450 MHz). The symmetric displacements in 
the velocity distribution, with a central hole and peaks on 
either side, indicates both deceleration and acceleration 
away from a center velocity. The observed velocity range 
is much less than expected for the normal BCF at this 
detuning. As mentioned previously, it was not possible 
to directly measure the acceleration or deceleration, but 
only their accumulated effects along the length of the 
slowing region. 

The center of the hole is at the center of the force pro- 
file, which for this measurement is Doppler shifted by 
800 MHz, or 493 7/fc. It cannot be caused by any of 
the individual bichromatic frequency components acting 
alone because each would have a large Doppler shift lo- 
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cated at either ±158 7/fc from the center of the velocity 
hole. We also confirmed that the effect is Dichromatic 
in nature by individually blocking one or more of the 
bichromatic frequency components, which resulted in a 
complete loss of the velocity shifts. 

It thus appears that a bichromatic force is present at 
2787, but that it is no longer a rectified force. This sug- 
gests that the sign of the force is changing sometime dur- 
ing the many stimulated cycles experienced by the atoms 
in between successive radiative decay events that reset 
the phase of the cycling. Indeed, we were able to simu- 
late the measured velocity profile quite well with a nu- 
merical Monte Carlo simulation that assumes a randomly 
directed force with a reduced magnitude that is typical 
of the BCF at 6I7 (25|. This supports a dephasing of the 
BCF, along with a reduction in magnitude caused cither 
by repeated sign changes or other effects. The previously 
observed vanishing of the force at 3687 suggests that the 
disruption of the BCF becomes even more severe at still 
larger detunings. In the next section we explore some 
possible explanations. 



C. Upper limits to the BCF 

We have examined a wide variety of physical and engi- 
neering limitations that could constrain the largest usable 
detuning for the BCF [25| , and here we describe the lim- 
iting factors most likely to arise in actual practice. One 
obvious issue is that the required laser power scales with 
S 2 , setting an engineering limitation. In terms of the sat- 
uration irradiance I a , the required irradiance from each 
beam (Eq. [5]) is 




A typical value of I s is a few mW/cm 2 , although for 
He* it is unusually small at 0.17 mW/cm 2 . Obtaining 
the necessary power is thus a significant challenge for cw 
experiments with detunings of several hundred times the 
natural linewidth 7. 

Another major consideration is the impact of imperfec- 
tions that cause cumulative dephasing of the stimulated 
cycling. The time-domain 7r-pulse model of Section III Al 
is particularly useful for estimating the size of these ef- 
fects. The self-cancellation of each pulse pair depends 
on precise matching of the left-hand and right-hand beat 
note "pulses" in Fig. [TJ If they are not symmetric, each 
successive pulse pair leaves the atom with a Bloch vector 
having a larger phase angle, and when the accumulated 
phase change reaches it, the direction of the force is re- 
versed. To make approximate estimates, we can replace 
the bichromatic beat notes with rectangular pulse trains, 
each with a repetition period T = tt/6 and pulse duration 
T/2, with their arrival times offset by T/4. If the reso- 
nant Rabi frequencies are Ox an( l ^2 for pulses arriving 
from the left and the right, the total accumulated phase 



error is approximately 

\ j* ^ + K(t)) 2 - ^ 2 + K(£ + T/4)) 2 ) dt. 

(9) 

Here the generalized Rabi frequencies incorporate possi- 
ble detunings uj[ and ui 2 of the bichromatic center fre- 
quencies from the atomic transition frequency, as in Eq. 
[21 The integration time t' is the interval between succes- 
sive resets of the stimulated cycling by radiative decay 
or collisional quenching to the ground state. It is given 
by t' = 4/7 for an ideal BCF cycle in which the atoms 
spend 1 /4 of their time in the excited state. 

The requirement for a well-defined BCF is that the 
dephasing in Eq. |H] is less than it. Because the Rabi 
frequency is proportional to the detuning, this becomes 
increasingly difficult to achieve at large detunings. Two 
possible sources of dephasing are: 

1. Intensity imbalance between the left-hand and right- 
hand bichromatic beams. This appears as a difference 
between f2i and Q 2 and leads to a first-order phase error. 
Setting the detunings to zero yields the condition Cli — 
il 2 < 2n/t'. For a nominal Rabi frequency of ft — 25 that 
yields a ir pulse in this approximation, the corresponding 
fractional imbalance is 

n 1 -n 2 _ c 7T l x -l 2 2-k 

— q— ~ *?• or ~r- ~ (10) 

using the proportionality of the irradiance / to f2 2 . In 
a badly disrupted BCF cycle the probability of atomic 
excitation will be close to the saturated statistical value 
of 1/2, so we set t' ~ 2/7 in Eq. [TO] This yields predic- 
tions generally consistent with experimental experience, 
although prior work on cesium suggests that the sensitiv- 
ity to beam imbalance is overestimated by about a factor 
of two [||. For He*, Eq. [10] predicts that a detuning of 
<5/27r=250 MHz requires irradiances balanced within 2%, 
and at 500 MHz, within 1%. In practice beam balancing 
beyond the 1% level is exceedingly difficult because of 
the imperfect spatial profile of the beams, accounting in 
part for the experimental observation of efficient slowing 
at 1857 (300 MHz) but not at 2787 or above. 

2. "Decelerative dephasing" due to rapidly chang- 
ing Doppler shifts as the atoms are slowed. A constant 
bichromatic force Ft, given by Eq. 0] will cause a time- 
dependent velocity change Av(t) — —Fbt/m, leading to 
detunings in Eq. [9] of +kAv(t) and -kAv(t + T/4). In 
Ref. [25[ a large effect from these detunings is estimated 
because the subtraction of the two terms in Eq. [S] is 
not fully taken into account. However, after this sub- 
traction only a very small quadratic term depending on 
T/4 remains, and no significant effect is predicted even 
for helium with its very small mass. Appreciable dephas- 
ing would occur only if there were also an overall shift 
ujq of the bichromatic center frequency from resonance. 
Then the detunings would become uj' x = ujq + kv(t) and 
uj' 2 = — kv(t + T/4), yielding a term proportional to 5. 
Asymmetric detunings ikv due to Doppler shifts, even 
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if large, do not cause this problem because they retain 
the mirror-image balance of pulses from the left and from 
the right. It is interesting to ask whether this cancella- 
tion is equally complete in the actual BCF configuration 
with its partially overlapping beat-note pulses, but ma- 
jor modifications to our present numerical modeling will 
be necessary to obtain a definitive answer. 

A different but somewhat related problem is the effect 
of rf phase errors between the right-hand and left-hand 
BCF beams that disrupt the nominal ir/2 phase shift be- 
tween successive beat notes. Phase errors can arise from 
external factors such as phase noise in the electronics, but 
also from the inherent variation of the rf phase with po- 
sition along the atomic beam axis. Numerical modeling 
of the OBEs is necessary to analyze the effects, because 
the quantitative dependence on the rf phase is not in- 
cluded in the 7r-pulse model. A full analysis is given in 
Section 3.5.2 of Ref. [25| . where it is shown that the 
BCF has a full width at half-maximum of about 9° in 
phase at 300 MHz detuning, or 5° at 600 MHz. At 300 
MHz the spatial length of an rf beat note is 50 cm, lead- 
ing to an unavoidable rf phase deviation of 14.4° across 
a slowing region 4 cm in length. For larger detunings 
this phase variation grows proportionally even while the 
phase sensitivity of the BCF grows worse, so little benefit 
is expected beyond this point. While this problem does 
not actually reverse the BCF deceleration, it reduces the 
force enough to make a single-stage slower for metastablc 
helium physically impossible. 

A final consideration is that at large irradiances, 
atomic energy levels can be modified by ac Stark shifts 
due to laser-induced coupling to distant states. For he- 
lium this is not a significant concern, as is shown in de- 
tail in Section 3.2 of Ref. (25|. However, in other atomic 
and molecular systems this could become a critical issue, 
since it not only leads to symmetric shifts uj[ and lo' 2 in 
Eq. [51 but also to rapid variations with time and position 
due to the overlapping traveling BCF beat notes. A brief 
discussion in the context of molecular slowing is given in 
Ref. 

Although we could not obtain sufficient experimental 
information at large detunings to definitively establish 
the cause of the anomalous result at 8 = 2787, it is clear 
that both intensity imbalance dephasing and rf errors are 
predicted to cause major problems at this detuning. To 
achieve velocity changes much larger than those observed 
at 5 = 1857, a different approach is needed. 



IV. CHIRPED BCF 

The chirped BCF decelerator accomplishes the same 
goal of slowing metastable helium atoms using bichro- 
matic forces without requiring the large bichromatic de- 
tunings ±S of the static slower described in Sec. IHI Al 
Instead we substitute a BCF with a relatively small de- 
tuning and thus a limited velocity range, in which the 
center frequency of each bichromatic beam pair is dy- 




FIG. 5. (color online). Three-dimensional BCF surface used 
in the Monte Carlo chirped BCF models. The beat note phase 
is assumed constant over the interaction region so that the 
force is a function only of the atomic velocity and Rabi fre- 
quency. The surface is interpolated from numerical solutions 
to the OBEs. The force is zero at Rabi frequencies less than 
0.45 and absolute velocities greater than IOO7/A;. 

namically adjusted to stay resonant with the He* atoms 
while they are slowed. This is done by linearly chirp- 
ing [29, 30] both laser frequencies to follow the changing 
Doppler shift. This greatly reduces the laser power re- 
quirements and allows the slower to operate well within 
the range of detunings for which the BCF is fully effective 
and reliable. 



A. Modeling 

To support our experiments, we have extended the nu- 
merical ODE calculations of Sec. Ill CI by using the force 
profiles to create a Monte Carlo model of chirped BCF 
slowing. This allows us to predict the performance as 
as a function of the chirp magnitude and rate. To ac- 
complish this, a two-parameter force function F(v, fi) is 
needed because both the atomic velocity and the Rabi 
frequency evolve slowly with time. The required profile 
is interpolated from a series of BCF profiles calculated 
for Rabi frequencies throughout the range from zero to 
y/3/25. The resulting BCF surface, depicted in Fig. [SJ 
shows a pair of prominent peaks at SI w 1.0 15 in addition 
to the primary maximum at SI = \j2>j1 S. 

Initial atomic velocities at the start of the Monte Carlo 
simulation are assigned at random using a probability 
distribution function derived from an empirical velocity 
distribution. The atoms are also assigned random initial 
time offsets between t = and —40 fxs to account for the 
finite atomic pulse duration, where t=0 corresponds to 
the beginning of the chopper aperture transit across the 
atomic beam path. The progress of an individual atom 
proceeds in small time steps and the change in velocity 
due to interaction with the bichromatic beams is evalu- 
ated at each step. The Monte Carlo model also accounts 
for the finite length of the chopped atomic beam pulse, 
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the gaussian bichromatic beam profile, the relative sizes 
of the atomic beam and bichromatic beam waists, and 
the transit time of the individual atoms into and through 
the bichromatic beams. A fuller account is given in Ref. 

Once all the atoms have passed through the interaction 
region, the time of flight for each atom is evaluated and 
the distributions of both the TOF and the final velocity 
are calculated to simulate experimental results. Our ini- 
tial results based on the full magnitude of the optimal 
BCF predicted that linear chirp ramps with durations of 
10 — 20 /j,s would work well, with little change in the size 
and shape of the slowed atom peak for chirps with mag- 
nitudes as large as 600 MHz. However, as we discuss in 
Sec. lIVBl in our experiments the optimal chirp durations 
were found to be in the range of 30 — 50 /is, with a clear 
dependence on the chirp magnitude. This was easily ac- 
counted for by reducing the effective BCF magnitude by 
a factor of two to approximately account for experimen- 
tal imperfections, at which point the Monte Carlo model 
predictions (shown in Fig. [5]) match the experimental re- 
sults quite well. 



B. Experiment 

The experimental layout for the chirped BCF exper- 
iments is shown in Fig. |6l highlighting differences from 
the large-detuning layout in Fig. [3] The chirped slower 
uses two Toptica DL100 diode lasers, one for the co- 
propagating +kv beam and the other for the counter- 
propagating — kv beam. This allows independent control 
of the co- and counterpropagating beam (±kv) Doppler 
shifts, a necessity for Doppler chirping. A weak sample 
beam from the —kv laser is frequency shifted by +kv 
using an AOM and sent to a saturated absorption spec- 
trometer where it is locked to the 2 3 Si — > 2 3 P2 transi- 
tion in helium. The other laser is stabilized at +kv by 
an offset locking scheme, in which the heterodyne beat 
note between the two lasers is measured on a photodi- 
ode, and locked by a microcontroller based phase locked 
loop (PLL) circuit to an offset of +2kv, typically cor- 
responding to about 1.6 GHz. The lasers are initially 
locked (without a chirp in progress) at the Doppler shift 
of a velocity group v(t ) near the center of the atomic 
velocity distribution, tog = uj a ± kv(to), and then chirped 
in opposite directions to follow the atomic Doppler shift 
as a function of time. Frequency modulation of the two 
lasers is accomplished using the a manufacturer-supplied 
model DL-MOD interface in one laser and a homemade 
copy in the other. The modulation is produced by an 
FET connected in parallel with the laser diode. A volt- 
age ramp applied to the FET gate causes a portion of the 
laser diode current to be diverted to ground, changing the 
laser output frequency. 

The bichromatic frequencies are generated for each 
laser separately using single-passed AOMs driven at an 
rf frequency 26. For each laser, the zero- and first-order 
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FIG. 6. (color online). Chirped BCF slower experiment block 
diagram showing key optical and electronic components. The 
configuration of the slower beams and atomic beam are the 
same as shown in Fig. [3] 



components are recombined using a 50% beam splitter 
to form the bichromatic beam. The beat note phase 
is controlled by locking the rf phase between a pair of 
homemade frequency synthesizer boards that drive the 
two AOMs. The microcontroller-based synthesizers share 
a common 10 MHz clock, and a small offset current is 
added to one of the PLL charge pump outputs to control 
the phase shift as described in Refs. [3ll [32j. 

We found that the rf phase was much less stable than 
expected because of microphonic motion of the lasers rel- 
ative to the AOMs, which causes phase variations be- 
cause of the short acoustic wavelength of the sound waves 
in the modulators. A feedback loop was added to provide 
additional phase stabilization by monitoring the relative 
rf phase of the optical beat notes using an Analog De- 
vices AD8032 phase detector. By using this phase mea- 
surement as an error signal in the PLL phase lock circuit, 
the RMS phase jitter was reduced to approximately 4°. 
However, frequent large phase excursions of up to ±28° 
could not be corrected and reduced the effective magni- 
tude of the BCF. A better solution might be to utilize a 
single AOM for both lasers, and this change is planned 
for future work. 

Our initial experiments were constrained to a very 
modest BCF detuning of 747 because of the limited op- 
tical power available from the DL100 diode lasers, which 
provided about 40 mW in each bichromatic beam pair 
after accounting for losses in the optics. For the same 
reason, the beams were focused to a top-hat radius of 440 
/jm that was somewhat smaller than the atomic beam, 
limiting the fraction of atoms that could be slowed. This 
situation could be greatly improved by the addition of a 
medium-power optical amplifier, such as a tapered laser 
amplifier diode. 
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FIG. 7. (color online). Experimentally observed velocity 
profiles for f requency-chirped BCF slowing with S = 747, 
Q r = a/3/2 5, and (f) — n/2. The four panels show chirp 
ramps of by MHz, 100 MHz in 30 [is, 200 MHz in 40 /is, 
and 300 MHz in 50 /is. 



C. Results and analysis 

The measured He* velocity distributions for chirped 
BCF slowing using a detuning of 747 and chirp mag- 
nitudes up to 300 MHz are shown in Fig. [7] The re- 
sults show the predicted increase in slowing with chirp 
magnitude, indicating that the chirped BCF profile re- 
mains resonant with a fraction of the atoms while they 
are slowed. 

An analysis of these chirp results is shown in Table UTJ 
in which the measurements of slowing and velocity range 
are defined in Sec. lIIIBl A detuning of only S — 747 with 
a 300 MHz frequency chirp provides a measured slowing 
of 210 7/fc or 370 m/s, the largest amount of slowing 
reported in any BCF experiment to date. This is more 
than twice the slowing measured for a static detuning 
of 5 = 1857 as reported in Table Q] The laser irradiance 
requirements are lower by nearly a factor of ten than 
what would be required for a static slower with the same 
velocity range, even if the problems outlined in Sec. IIII CI 
could be overcome. 

The optimal experimental chirp ramp durations are 
found to increase from about 30 to 50 fj,s as the chirp 
magnitude increases from 100-300 MHz. As mentioned 
previously, the Monte Carlo model described in Sec. lIV Al 
yields results in good agreement if the bichromatic force 
magnitude is reduced by a factor of two from its ideal 
value. This reduction could easily be caused by experi- 
mental imperfections such as rf phase jitter and imperfect 
gaussian beam profiles. 

The full velocity profiles predicted by these model cal- 
culations are shown in Fig. [5J Again the agreement 
with experiment is good. The most probable velocity of 



TABLE II. Summary of experimental BCF chirp results for 
S — 747. All experiments used the optimal BCF parameters 
Q r ~ \/3/2 5 and 4> = t/2. Slowing is reported both in units 
of 7/fc = 1.76 m/s and 5/k = 130 m/s to facilitate comparison 
with static BCF deceleration. 
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FIG. 8. (color online). Chirped BCF Monte Carlo model out- 
put after reducing BCF magnitude, using same experimental 
parameters used to obtain the results in Fig. [7] 



the chirped group of atoms experimentally matches the 
model prediction within 20 j/k, although the fraction 
of slowed atoms is relatively small in this initial exper- 
imental configuration. There was no slowing observed 
when the chirp was increased to 400 MHz, indicating 
badly sub-optimal BCF conditions. The model predicts 
this behavior in two scenarios — a smaller BCF magni- 
tude than expected, or periodic, temporary loss of the 
BCF as would be caused by occasional large phase er- 
rors. Both of these factors were potentially present in 
the experimental configuration for this initial test. 

The small fraction of slowed atoms was caused not only 
by limited beam overlap, but by practical difficulties in 
using a chopped atomic beam to measure velocity distri- 
butions with a single frequency-chirp ramp. The atomic 
beam pulse duration, about 40 /is at the chopper, ex- 
pands to more than 120 [is by the time the atoms reach 
the interaction region. This is longer than the chirp ramp 
duration, so atoms on the leading edge and trailing edges 
of the pulse fall outside the chirped BCF profile and are 
not slowed. This effect is more pronounced with faster 
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frequency ramps and results in a much reduced number 
of slowed atoms. In a fully realized design, rapidly re- 
peating chirp sequences can be used in conjunction with 
a cw atomic beam to assure that most of the atoms in a 
continuous beam experience the full chirped BCF force, 
as discussed in Sec. |V] 



V. CONSIDERATIONS FOR PRACTICAL 
ATOMIC DECELERATION FOR MOT LOADING 

While initial attempts to slow He* atoms using the 
chirped BCF scheme described in Sec. HVBl were limited 
to Av < 370 m/s, we were able to identify some easily 
remedied shortcomings in the experimental layout and 
verify the accuracy of the Monte Carlo model discussed 
in Sec. HVAl 

Even without major changes, simply scaling the 
300 MHz chirp at 8 = 747 up to a detuning of approxi- 
mately 1237 and better matching the bichromatic beam 
to atomic beam diameters will permit a frequency chirp 
magnitude of 600 MHz. This chirp, assuming an initial 
Doppler shift of 800 MHz, is the minimum required to 
slow a reasonable number of atoms to a final velocity of 

50 j/k. Reducing the initial Doppler shift would per- 
mit a smaller frequency chirp magnitude at the expense 
of the number of slowed atoms. 

The brightness of such a chopped atomic beam slower 
would be very low due to the atom beam duty cycle of 
only 1%. To create a slower capable of replacing a Zee- 
man slower, the chirp configuration must be adapted to 
use a continuous or near-continuous atomic beam. This 
requires that the frequency chirping cycle repeats con- 
tinuously and as quickly as possible. Using a chirp ramp 
duration of 10 fjs with an additional 5 /is to reset the 
laser frequency gives a cycle period of 15 /xs. However, 
because the 5 fis is much less than the atomic transit 
time through the slower, it will not appreciably interfere 
with the cw operation of the slower. The laser frequency 
locking must occasionally be re-established, but it suf- 
fices to do this once every few milliseconds. For these 
estimates, we will assume the laser locks require 500 fis 
to reset, which must occur every 2 ms. This results in a 
fairly realistic estimate of an 80% duty cycle. 

We estimate the output flux of such a slower to facili- 
tate comparisons to other atomic beam slowing methods. 
Assuming conservatively that 10% of the metastable he- 
lium atoms are subject to slowing after taking into ac- 
count the beam overlaps, BCF velocity range, and initial 
Doppler shift, we can predict the intensity of slow atoms 
using information from the model simulation shown in 
Fig. [5] From Sec. IIII Al the metastable source bright- 
ness is 3 x 10 13 He* s _1 Sr -1 . Adjusting for the duty 
cycle and slower efficiency, the BCF slowed atom bright- 
ness will therefore be 0.8 x 0.1(3 x 10 13 He* s^ 1 Sr" 1 ) 
2.-1 x 10 12 He* s _1 Sr -1 . Factoring in the measured ac- 
ceptance angle of the slower and its estimated exiting 
beam diameter of 500 /im, the estimated slow atom flux 
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FIG. 9. (color online). Predicted velocity distribution of a 
fully realized chirped BCF slower acting on a cw atomic beam. 
The black dashed line is the velocity distribution of the atomic 
beam source, and the red solid line is the predicted distribu- 
tion with the chirped BCF present. Notice the large slow 
atom peak at ~ 5O7/A; and relatively uniform depletion of a 
large range of atomic velocities. 



is (2.4 x 10 12 He* s" 1 Sr~ 1 )(1.5 x 10~ 5 Sr)/0.2 mm 2 ~ 
2 x 10 s He* s^mm -2 . 

This can be compared with the 1083 nm Zeeman slower 
used by Vassen group, which has a MOT loading time of 
0.5 s [33l439| . Their scheme begins with a dc discharge 
He* source somewhat similar to ours, but with a higher 
brightness of 3 x 10 14 He* s _1 Sr -1 . A two-part Zee- 
man slower [33ll37| is used in conjunction with an optical 
molasses pre-collimator that increases the angular accep- 
tance of the slower by a factor of ~ 30 [36[ . The output 
flux of slow He* atoms is 2 x 10 9 He* s" 1 mm" 2 W\. Tak- 
ing into account the brighter initial metastable helium 
source, this is the same flux that would be predicted for 
a chirped BCF decelerator without a pre-collimator. The 
main reasons for the improved performance of the BCF 
slower are the greatly reduced length and the absence of 
transverse heating effects common in long Zeeman show- 
ers. With the addition of a pre-collimator the output flux 
could be further increased by a factor of 2-50, depending 
on the method used [la [39|, |40j ■ 

Development of an improved chirped slower able to 
slow a cw atomic beam to rest is currently underway in 
our laboratory. 

VI. SUMMARY 

Our investigation of the bichromatic force at large 
detunings shows that a static two-stage slower for a 
metastable helium beam is probably not feasible. We 
have added to prior experimental evidence for an effec- 
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tive upper limit to the BCF at a detuning of about 2507 
for He*, probably due to accumulated small dephasings of 
the stimulated Rabi cycling in between "resets" by radia- 
tive decay. In addition, it is difficult to reliably produce 
the required laser irradiance at larger dctunings because 
it increases quadratically. 

To circumvent these limits a novel BCF slower was 
developed, using a small detuning while chirping the 
Doppler shift of the force center to maintain resonance 
with the decelerating atoms. In a prototype experi- 
ment we have achieved slowing of metastable helium by 
210 7/fc (370 m/s) with a detuning of only 747. This 
is more than twice the largest BCF slowing previously 
attained, and represents a savings in laser irradiance of 



nearly a factor of 10. Scaling up the detuning to 1237, 
increasing the bichromatic beam waists, and rapidly re- 
peating the chirp will allow the realization of a He* slower 
capable of slowing atoms to MOT capture velocities. The 
chirped BCF slower is predicted to have a brightness at 
least comparable to current Zeeman slowers, but with a 
much shorter length. 
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